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FOREWORD 


This  report  was  prepared  by  Textron's  Bell  Aerosystems  Company,  Buffalo, 
New  York,  under  Contract  AF33(657)-8936.  The  work  was  administered  under  the 
direction  of  the  Flight  Dynamics  Laboratory,  Aeronautical  Systems  Division,  by 
Mr.  G.  E.  Maddux,  Project  Engineer.  The  work  was  performed  by  the  Structures 
Section  of  the  Aerospace  Engineering  Department,  Bell  Aerosystems  Company, 
in  the  period  of  15  June  1962  to  31  July  1963.  Mr.  Richard  H.  Gallagher  was 
Technical  Director  of  the  study. 

The  authors  wish  to  acknowledge  the  contributions  to  this  work  of  P.P. 
Bijlaard,  Professor  of  Engineering  Mechanics  at  Cornell  University  and  consultant 
to  the  Bell  Aerosystems  Company.  Prof.  Bijlaard  formulated  the  techniques  em¬ 
ployed  to  obtain  solutions  for  the  instability  of  flat  sandwich  panels  and  for  the 
thermal  stresses  and  instability  of  sandwich  cylinders  supported  by  rigid  bulkheads. 
Major  contributions  to  this  work  were  also  made  by  Mrs.  Leona  Barback,  and 
Miss  Penelope  Miller  of  the  Bell  Electronic  Data  Processing  Department,  who 
coded  the  computer  programs  used  in  development  of  the  design  data  presented 
herein. 
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ABSTRACT 


Structural  design  data  are  presented  for  various  applied  load  and  thermal 
gradient  conditions  for  flat  rectangular  sandwich  panels  and  sandwich  cylinders. 

The  flat  panel  solutions  pertain  to  instability  under  nonuniform  stress,  and  also  to 
the  stresses  and  displacements  resulting  from  normal  pressures  and  temperature 
gradients  across  the  panel  thickness  in  the  presence  of  uniform  midplane  com¬ 
pression. 

Sandwich  cylinder  design  data  are  given  for  buckling  under  nonuniform  circum¬ 
ferential  and  axial  stress,  respectively,  and  for  stresses  due  to  radial  or  axial 
temperature  gradients.  The  range  of  stiffness  parameters  extend,  at  one  limit,  to 
the  isotropic  forms  of  construction. 


PUBLICATION  REVIEW 
This  report  has  been  reviewed  and  is  approved. 


Colonel,  USAF 


Chief,  Structures  Division 


AF  Flight  Dynamics  Laboratory 
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CHAPTER  I 
INTRODUCTION 


The  development  of  thermal  stress  determination  techniques  for  supersonic 
transport  aircraft  structures  has  consisted  of  a  threefold  effort.  First,  a  survey  of 
the  literature  pertinent  to  thermal  stress  analysis  was  performed,  resulting  in  pub¬ 
lication  of  an  annotated  bibliography  of  the  subject  for  the  period  1955-1962  (Reference 
1);  a  bibliography  covering  thermal  stress  analysis  literature  prior  to  1955  had  been 
published  elsewhere.  The  second  portion  of  this  study  effort  has  been  concerned  with 
the  development  of  graphical  design  data  for  various  applied  load  and  thermal  gradient 
conditions  for  flat  rectangular  sandwich  panels  and  sandwich  cylinders.  This  work  is 
presented  in  the  subject  report.  Finally,  computer  programs  for  the  more  complex 
problems  encountered  in  beam,  plate,  and  cylindrical  constructions  were  coded  with 
the  purpose  of  making  them  available  to  eligible  recipients  in  the  airframe  industry. 

A  description  of  these  programs  is  given  in  Reference  2. 

Sandwich  panels  and  cylinders  provide  attractive  constructional  forms  for 
large  high-performance  vehicles  and  will  undoubtedly  be  given  consideration  during 
the  developmental  phase  of  the  SST  airframe.  They  are  extremely  efficient  from  a 
weight  standpoint  and  can  be  fabricated  using  materials  suitable  for  an  elevated 
temperature  environment  exceeding  1000°  F.  Sufficient  data  have  been  accumulated  to 
permit  their  design  for  conventional  loadings  and  uniform  temperature  conditions 
(c.f.,  Reference  3)  but  design  data  is  lacking  for  nonuniform  temperature  conditions. 

It  is  the  intent  of  this  report  to  provide  such  data.  In  addition,  it  is  to  be  noted  that 
the  design  curves  to  be  presented  extend,  in  one  limit,  to  the  conventional  isotropic 
thin  plate  and  cylinder  constructions. 

The  design  curves  are  presented  in  the  chapters  of  this  report  with  a  dis¬ 
cussion  of  their  basis,  use,  and  limitations  and  with  illustrative  examples,  but  without 
a  detailed  description  of  the  related  formulations  and  solution  techniques.  The  latter 
are  given  in  a  series  of  Appendixes. 

Solutions  to  flat  sandwich  panel  problems  are  found  in  Chapters  II  and  III.  One 
of  the  most  important  problems  in  flat  plate  analysis  is  the  prediction  of  elastic 
instability  in  the  presence  of  nonuniform  stress  produced  by  temperature  gradients. 
HoffW,  Klosner  and  Forray  and  van  derNeut  ^ ,  among  others,  have  presented 
procedures  and  isolated  results  for  this  problem  as  it  pertains  to  thin  isotropic  plates; 
there  are  no  known  solutions  for  sandwich  construction,  In  Chapter  II  of  this  report, 
curves  are  presented  which  allow  the  prediction  of  the  elastic  instability  conditions 
for  long  honeycomb  sandwich  panels  with  equal  thickness  isotropic  skins  under  non- 
uniform  longitudinal  stress,  for  both  fixed  and  simply  supported  longitudinal  edges. 

Manuscript  released  by  authors  in  July,  1963  for  publication  as  an  ASD 
Technical  Documentary  Report. 
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Chapter  III  is  coneerned  with  the  stresses  and  displacements  at  the  eenter  of 
simply-supported  rectangular  sandwich  panels  subjected  to  temperature  gradients 
through  the  panel  thiekness  and  to  pressure  loadings.  The  presence  of  fixed  uniaxial 
midplane  forees  is  also  taken  into  aeeount.  The  graphs  eover  a  wide  range  of  aspect 
ratios,  stiffnesses,  and  midplane  loadings.  A  limited  number  of  results  for  conditions 
of  fixed  support  along  two  edges  is  also  presented.  Techniques  for  solving  this  type 
of  problem  have  been  published  by  Bijlaard  (^)  and  Ebeioglu  (8).  Reference  7  con¬ 
siders  only  the  problem  of  temperature  gradients  and  in  neither  Reference  7  nor  in 
Reference  8  are  design  eurves  presented.  The  formulation  of  the  governing  dif¬ 
ferential  equations  for  this  problem,  as  described  in  Appendix  B,  has  not  been 
formally  presented  elsewhere. 

Chapters  IV  through  VII  deal  with  various  problems  in  the  analysis  of  heated 
sandwieh  cylinders.  Curves  for  the  prediction  of  the  elastie  instability  of  sandwieh 
cylinders  subjected  to  cireumferentially  varying  axial  stress  are  found  in  Chapter 
IV.  Conditions  which  produee  this  type  of  stress  distribution  inelude  both  nonuniform 
eireumferential  temperature  gradients  and  also  combinations  of  bending  and  axial 
compression.  Graphical  means  of  determining  the  stresses  due  to  temperature  for 
this  problem  have  not  been  considered  since  the  applicable  procedures  involve  the 
use  of  either  thermoelastie  beam  theory  or  the  diserete  element  techniques  of 
matrix  structural  analysis.  In  neither  ease  are  graphical  solutions  feasible. 

Published  references  have  not  as  yet  considered  the  above-cited  instability 
problem.  Bijlaard  and  Gallagher  (  ® )  and  Abir  and  Nardo  studied  the  ease  of 
the  isotropie  cylinder  under  circumferentially  varying  axial  stress.  Both  references 
concluded  that  the  small  defleetion  theory  maximum  stress  amplitude  for  elastie 
instability  under  varying  stress  is  effectively  equal  to  that  for  uniform  compression. 
This  eonelusion  is  found  to  be  essentially  correet  in  the  present  ease  of  sandwich 
cylinders.  Comparisons  are  made  with  the  results  of  past  studies  of  the  stability 
of  sandwieh  cylinders  under  uniform  compression  (References  11-13). 

Chapters  V  and  VI  treat  the  stress  and  instability  analysis  of  heated  cylinders 
supported  by  unheated  rigid  bulkheads.  Fixed  and  simple  support  conditions  are  in- 
eluded.  Graphical  methods  for  determining  the  maximum  stresses  under  such 
conditions  appear  in  Chapter  V.  This  class  of  problem,  for  isotropic  cvlinders  has 
been  previously  dealt  with  by  Hoff^\  Przemienieeki  and  Johns  Based  on 

the  results  of  Chapter  V,  a  series  of  elastic  instability  analyses  were  performed  and 
from  these  latter  results  a  set  of  graohical  representations  of  the  data  were  drawn. 
These  appear  in  Chapter  VI.  Results  for  the  limiting  case  of  the  isotropic  cylinder 
are  compared  with  solutions  given  by  Hoff  (^)f  Johns  and  Anderson  (1®. 

The  final  design  data  given  in  this  report  pertain  to  the  stresses  in  a  sandwich 
cylinder  due  to  radial  temperature  gradients.  Means  for  computing  the  stresses  due 
to  either  internal  or  external  pressure  are  also  presented.  A  similar  problem,  but 
with  different  core  properties,  was  treated  by  Yao  Yao  proposes  a  numerical 
approach  to  analysis,  based  on  a  method  of  successive  approximations.  The  radial 
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temperature  gradient  problem  is  limited  to  stress  solutions  since  considerations  of 
general  instability  are  involved. 

As  noted  earlier,  all  theoretical  questions  associated  with  the  computation  of 
the  data  employed  in  construction  of  the  graphical  representations  are  examined  in 
detail  in  a  series  of  Appendixes  (Appendixes  A  through  E). 
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CHAPTER  II 

BUCKLING  OF  SANDWICH  PANELS  UNDER  NONUNIFORM  STRESS 


This  chapter  presents  graphical  representations  with  which  the  elastic 
instability  conditions  for  honeycomb  sandwich  panels  with  equal  thickness  isotropic 
skins  under  nonuniform  longitudinal  stress  can  be  predicted,  The  panels  are 
assumed  to  be  "long"  (i,c.,  aspect  ratio  effects  and  the  significance  of  the  trans¬ 
verse  support  conditions  are  excluded)  and  clamped  or  simply  supported  along 
their  longitudinal  edges.  The  conditions  of  analysis  are  illustrated  in  the  insert 
in  Figure  II-l. 

A  detailed  treatment  of  the  formulations  and  procedures  employed  in  the 
derivation  of  the  presented  curves  is  given  in  Appendix  A.  The  governing 
differential  equation  is  that  which  was  derived  by  Reissner  in  Reference  20;  the 
method  employed  for  its  solution  was  the  finite  difference  approximation  technique. 
A  special  purpose  computational  program  was  coded  and  utilized  in  the  development 
of  the  required  data. 


It  is  shown,  in  Appendix  A,  that  the  critical  stress  is  a  function  of  two  para¬ 
meters,  a  stiffness  parameter  (D  )  and  a  parameter  defining  the  stress  distribution. 
The  stiffness  parameter  is  given  *as 


D 

P 


(II- 1) 


where  b  is  the  panel  width  and  D  and  D  are  the  core  shear  stiffness  and  the  panel 
bending  stiffness,  respectively.  C*The  analyst  has  at  his  disposal  the  choice  of  a! 
number  of  formulas  for  Dq  and  Ds,  dependent  upon  the  details  of  sandwich  panel 
construction;  these  choices  are  developed  and  discussed  in  Reference  3,  Section  3.1. 
The  simplest  and  most  appropriate  formulas  for  these  parameters  are 


Dq  (hM)Gc 

n  _  Et  (h+t)2 

s  2(1-^  2) 


(II- 2) 

(II— G) 


where  h  is  the  core  depth,  G  is  the  core  modulus  of  rigidity,  t  is  the  thickness  of 
one  of  the  panel  faces,  and  ffand  jj.  are  the  modulus  of  elasticity  and  Poisson's 
ratio  of  the  face  material,  respectively.  These  designations  of  geometry,  stiffness, 
material  properties,  etc.  are  applicable  to  all  chapters  of  this  report. 


In  the  strictest  sense,  the  stiffness  paratmeter  Dp  can  range  from  nearly 
zero  to  infinity.  At  infinity  'he  case  of  the  isotropic  plate  with  infinite  shear  stiff¬ 
ness  prevails.  This  conditi  ;ii  is  nearly  achieved  at  much  lower  values  and  it  is 
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satisfaetory  to  eonsider  the  maximum  Dp  to  equal  10  .  At  the  lower  range  of  Dp  it 
is  found  that  an  entirely  different  mode  of  instability,  wrinkling,  is  encountered. 

This  mode  is  characterized  by  a  zero  wavelength  buekling  condition  due  to  an  initial 
assumption  that  the  individual  bending  stiffnesses  of  the  faees  eould  be  neglected.  It 
is  not  the  intention  of  the  present  analysis  to  deal  with  wrinkling.  The  curves  have 
consequently  been  terminated  at  ^x  =  1,  whieh  is  a  limit  imposed  by  wrinkling 

(h+tlGc 

associated  with  shear  instability  failure  of  the  eore.  Thus,  evaluations  of  the 
critical  stress  were  performed  for  a  range  of  Dp  values  from  10**  down  to  values 
where  this  shear  instability  limit  appeared. 

Only  the  simplest  forms  of  stress  distribution  are  of  general  interest  and 
therefore  suitable  for  inclusion  in  design  charts.  The  results  presented  herein  were 
developed  for  linear  and  trigonometric-shaped  stress  distributions.  Linear  or  nearly 
linear  distributions  exist  on  panels  in  wings  sustaining  overall  nonuniform  chordwise 
temperature  gradients  and  also  when  the  bending  stress  distribution  on  a  wing  cross- 
seetion  departs  from  the  uniformity  predicted  by  elementary  theory  (as  is  generally 
the  case).  Trigonometric-shaped  stress  distributions  are  good  approximations  to 
conditions  where  the  panel  longitudinal  edge  members  (e.g..  the  spar  caps)  provide 
significant  heat  sinks,  thereby  resulting  in  extreme,  loeal,  transverse  variations  of 
the  longitudinal  thermal  stress.  Similar  comments  as  to  applicability  could  be  made 
for  both  types  of  stress  distribution  for  fuselage  skin  panels. 

A  linear  variation  of  stress  which  is  equivalent  to  a  uniformly  distributed  load 
superimposed  upon  a  bending  moment  is  conveniently  described  by  the  ratio  (S)  of  the 
stress  levels  at  either  side  of  the  panel.  Using  the  convention  that  the  denominator 
is  always  the  larger  positive  (compressive)  value,  any  purely  eompressive  load  is 
covered  by  the  range  1  —  S  -  0.  This  range  has  been  extended  into  the  mixed 
compressive-tensile  region  to  the  value  S  =  -1.0,  i.e.  pure  bending. 

This  range  of  S  from  -*-1.0  to  -1.0  has  been  covered  in  intervals  of  0.25.  The 
selected  trigonometric  cases  are  simply  sin  y,  sin  2y,  eos  y,  cos  2y. 

Figures  II-l  to  II -4  present  the  results  in  graphical  form.  Figures  II-l  and 
II-2  are  for  the  simply  supported  case  while  Figures  11-3  and  11-4  relate  to  fixed 
support  of  the  longitudinal  edges.  It  was  intended  to  depict  in  Figure  II-l  the  results 
from  Reference  3  for  uniform  compression  (S  =  1.0),  but  due  to  the  excellent  level  of 
agreementitis  not  possible  to  differentiate  between  those  results  and  the  present 
solutions  for  S  =  1.0.  The  eurves  are  terminated  at  the  lower  end  by  the 
wrinkling  boundary  diseussed  previously. 

Figure  II-5  depicts  the  variation  of  buckling  stress  for  the  isotropic  plate  sub¬ 
ject  to  linear  stress  distributions  for  each  boundary  condition.  These  curves  have 
been  obtained  by  considering  large  values  of  Dp.  Here,  the  form  of  thcabeissa  has 
been  altered  slightly  to  eliminate  the  shear  stiffness  of  the  core  (D().  The  new 
abcissa  is  thus  taken  as  M  A*  cr  b2  t  where  tj  2t  is  the  thickness  if  the  isotropic 
plate. 
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Illustrative  Examples 


Example  (1)  Sandwich  Plate  with  simply  supported  edges 


Width  of  panel 
Core  Thickness 
Modulus  of  faces 
Face  Thickness 
Shear  Modulus  of  Core 


=  40  in. 

=  0.50  in, 

107  lb/in? 
0.036  in. 

=  2  x  104  lb/in.2 


The  panel  is  subjected  to  a  triangular  compressive  load  (S  =  0).  Utilizing  the 
above  data  and  formulas  3.12E  and  3.13B  from  Reference  3,  one  obtains  Ds  =  5.68  x 

D 

104  and  D  =  1.0  x  104.  From  this  paper,  D  =  <4  =  282.  From  figure  II-l,  for 

4  ^  r\ 


2 1 

S  =  0  at  Dp  =  282,  the  buekling  load  parameter  -jy—  <r  is  found  to  be  0.257.  Henee, 

0.257  D  q 

<J  =  - — ~C|  =  35700  lb/in2 

or  21 


Total  applied  load  -  2t  cr  —  =  51400 

cr  2 


lb. 


Example  (2)  Isotropic  Plate  with  simply  supported  edges 


Width  of  panel 
Modulus  of  panel 
Thickness 
Thus  I) 

s 


=  40  in. 

107  lb/in? 
0.388  in. 
5.343  x  104 


This  panel  is  also  subjected  to  a  triangular  compression  (S  =  0).  From  Figure 
tjb^ 

II-5  at  S  =  0.  cr  =  77.05 

Ds  er 

77.05  D  .  2 

Thus  <j  -  _ s  =  6625  Ib/in. 

cr"  ,  ,  o 


Here  again  total  applied  load  =  tj  51400  lb.  Thus,  the  two  panels 

carry  the  same  loading  but,  as  can  be  seen,  the  sandwich  panel  has  a  considerable 
advantage  over  the  isotropic  plate  on  a  weight  basis. 


The  few  published  solutions  for  budding  of  long  panels  under  nonuniform  load¬ 
ing  are  entirely  restricted  to  the  case  of  isotropic  plates.  Timoshenko  (21)  has 
considered  an  isotropic  plate  with  a  linearly  varying  stress.  The  parameter  used 
by  Timoshenko  is  the  factor  k  in  the  expression. 


a 

cr 


tt^Dc 

2 

b  t 
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Using  this  parameter  the  following  comparison  is  obtained 


s 

k 

Timoshenko 

Present  Report 

+  1.0 

4.00 

3.98 

0 

7.81 

7.79 

-0.5 

13.40 

13.35 

-1.0 

23.90 

23.63 

(4) 

Hoff  has  considered  a  trigonometric  distribution  of  load,  corresponding  to  the 
cos  2 y  case  of  the  present  paper.  The  value  obtained  by  Hoff,  k  =  7.67,  shows  good 
agreement  with  the  present  work  where  k  =  7.68. 


S 


Stiffness  Parameter  D 


a 


Buckling  Load  Parameter  FT  =  2t  <r  /D 

cr  q 


Figure  II-l.  Buckling  Load  of  a  Long  Sandwich  Panel  Subjected  to  Linearly 
Varying  Edge  Stress  (Long  Edges  Simply  Supported) 


Figun 
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I 


I 

Figure  II-5.  Buckling  Stress  for  Isotropic  Pla 
Varying  Compressive  Load: 
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CHAPTER  III 

PROCEDURE  FOR  DETERMINING  THE  STRESSES  AND 
DISPLACEMENTS  FOR  A  RECTANGULAR  SANDWICH  PANEL 


Design  curves  are  presented  in  this  chapter  which  permit  the  determination  of 
the  central  deflection  and  maximum  bending  moments  in  a  rectangular  sandwich 
plate  subjected  to  a  uniform  normal  pressure  and  uniform  temperature  gradient 
through  its  thickness  combined  with  an  unaxial  in-plane  compression  Nx.  All  edges 
of  the  panel  are  simply  supported.  In  addition  a  limited  number  of  curves  is 
presented  for  the  case  of  two  opposing  edges  fully  built-in  and  two  edges  simply 
supported.  The  sandwich  panel  consists  of  two  isotropic  faces  of  equal  thickness 
t(Fig.  III-l)  and  a  core  of  thickness  h.  The  core  has  the  usual  properties  of  zero 
in-plane  stiffness  and  uniform  transverse  shear  stiffnesses  with  a  shear  modulus 
Gc. 


The  type  of  loading  combination  considered  here,  consisting  of  mechanical 
and  thermal  loading,  can  occur,  for  example,  in  a  wing  panel  of  a  high  speed  air¬ 
craft  subjected  to  kinetic  heating  effects  and  cooled  by  the  presence  of  fuel  stored 
in  integral  tanks.  The  values  of  the  mechanical  and  thermal  loadings  are  assumed 
to  have  been  determined  from  a  general  structural  analysis. 

The  governing  small  deflection  theory  differential  equations  for  a  sandwich 
panel  subjected  to  thermal  as  well  as  mechanical  loading  have  been  presented  by 
Bijlaard(?)  and  Ebcioglu(S)  but  in  neither  paper  design  curves  are  presented.  In 
Reference  22  a  new  formulation  of  the  differential  equations  is  developed  and  the 
equations  are  solved  for  a  number  of  boundary  conditions.  An  outline  of  the  method 
of  solution  is  given  in  Appendix  B.  The  solution  which  involves  the  use  of  infinite 
Fourier  series  has  been  programmed  in  FORTRAN  for  evaluation  on  an  IBM  7090 
computer. 

As  in  Chapter  II  the  stiffness  of  the  sandwich  panel  is  described  by 

single  pai'ameter  2 

D  b  D 


s 


In  addition,  however,  since  panels  of  finite  proportions  are  now  being  considered, 
the  aspect  ratio 

A 

must  also  be  introduced.  Although  both  Dp/  and  A  may  both  vary  over  very  large 
ranges  theoretically, practical  considerations  provide  limitations,  so  that  only 
values  of  Dj^  between  1  and  100  and  A  between  0.4  and  2.5  have  been  considered. 
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Since  there  are  three  independent  parameters,  stiffness,  aspect  ratio  and  mid-plane 
loading,  it  is  not  possible  to  present  all  variations  simultaneously  on  one  graph. 
Thus  Figures  III-l  to  III-4  present  the  variation  of  the  maximum  deflection  per 
unit  temperature  gradient  for  four  values  of  the  mid-plane  loading  with  all  edges 
simply  supported.  Similarly  Figures  III-5  to  III-8  show  curves  of  the  maximum 
deflection  per  unit  normal  pressure  for  the  same  four  mid-plane  loadings  with  all 
edges  simply  supported.  Figures  III-9  to  III- 16  depict  the  central  bending  moments 
per  unit  temperature  difference  or  unit  pressure  for  the  four  mid-plane  loadings 
with  all  edges  simply  supported. 

In  Figures  III- 1 7  and  III- 18  the  central  deflections  per  unit  temperature 
gradient  and  per  unit  normal  pressure  are  plotted  against  the  mid-plane  loading  for 
a  panel  having  a  stiffness  Dp'  =  100  (i.e.  an  equivalent  isotropic  plate).  For  values 
of  the  aspect  ratio  up  to  unity  the  buckled  shape  of  such  a  plate  consists  of  a  single 
half-wave  in  each  direction  and  thus  buckling  is  represented  by  the  maximum  de¬ 
flection  increasing  smoothly  to  infinity.  On  the  other  hand,  for  a/b  2  1.5  the  buckled 
shape  has  two  or  more  half  waves  in  the  longitudinal  direction  and  therefore  does 
not  conform  to  the  deflection  form  generated  by  the  thermal  or  normal  loading.  The 
curves  for  a/b  =  1.5,  2.0  and  2.5  have  been  terminated  at  the  true  values  of  the 
critical  stress  as  given  by  Timoshenko  for  an  isotropic  plate.  At  these  values  the 
buckling  will  be  characterized  by  a  sudden  change  in  deflected  shape  of  the  panel. 
Similar  curves  may  be  generated  by  suitable  cross  plotting  for  other  values  of  Dp'. 

Figures  III-19  to  III-22  show  central  deflections  and  bending  moments  similar 
to  Figures  III-l,  III- 5,  III- 9  and  III- 13,  but  with  two  opposing  edges  fully  built-in 
and  two  edges  simply  supported.  The  bending  moments  at  the  center  of  the  built-in 
edge  are  given  in  Figure  III-23. 

Illustrative  Examples 

(1)  A  rectangular  sandwich  panel  20  in,  by  30  in,  is  simply  supported  along 
each  of  its  edges  and  is  subjected  to  a  uniform  normal  pressure  of  1  psi.  The 
dimensions  of  the  sandwich  are  as  follows: 

Core  thickness  =  0,50  in.  E  =  10.5  x  106  psi 

Face  thickness  =  0.020  in.  Gc  =  6.65  x  104  psi 

Using  the  above  data  D  =  3.15  x  104  and  D  '  =  44.5.  By  interpolation  on  Figure  III- 
5  for  Dp'  =  44.5  and  AS  =  1.5 


0.77 

w  =  0.0401  in. 
max 


4^ 

w  IT  D„ 
max  s 


pb 

Hence  maximum  deflection 
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In  Reference  23  a  sandwich  having  the  above  dimensions  has  been  sub¬ 
jected  to  normal  pressure.  The  maximum  deflection  measured  on  this  specimen 
was  0.0384  in.,  compared  with  a  predicted  maximum  deflection,  using  the  theory 
of  Reference  23,  of  0.0446  in. 


(2)  A  rectangular  sandwich  panel  24  in,  by  36  in,  is  simply  supported  along 
each  edge  and  is  subjected  to  a  uniaxial  compression  of  1500  lb/in.  The  panel  also 
carries  a  normal  pressure  of  10  psi  acting  downwards.  The  lower  face  is  heated  so 
that  there  is  a  constant  temperature  difference  of  200 °F  between  the  lower  and 
upper  faces.  It  is  required  to  find  the  maximum  deflection  at  the  center  of  the  panel 
and  also  the  maximum  bending  moments. 

Core  thickness  =  0.97  in.  E  —  10.5  x  10®  psi 

Face  thickness  =  0.03  in.  G  =  3  x  10^  psi 

ct  =  1.2  x  10-5  /  °F 

Using  the  above  dimensions  Ds  =  1.73  x  105,Dq  =  3  x  104,  A  =  1.5 

N 

x 

Thus  D  '  =  10,0  and  rj  -  ~ —  =  -0.05 
P  x  Dq 

From  Figure  III-2  for  these  values  of  Dp'  and 
_  2 

w  7T  h 
max 

— -  ■ 1135 

b  (1+  /i  )  T 


w  =  0.2067  in.  (downwards) 
max 

T 


From  Figure  III- 6 


a 

w  ir  D 
maxp  £ 

Pb4 


=  0.97 


Wmax  =  °«1910  (downwards) 


Total  maximum  deflection  =  0.2067  +  0.1910  =  0.3977  in. 


In  a  similar  fashion  using  Figures  III-10  and  III- 14  the  maximum  bending 
moments  are  found  to  be 

M.\  Mxp  <-  Mxt  =  294.14  -  250.45  =  43.69  lb.  in. 

M  Myp  +  M  t  =  535.75  -  16.19  =  519.56  lb.  in. 
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These  bending  moments  yield  direct  stresses  in  the  facings  crx  -  ±1262  psi  and 
J1  =  ±15015  psi.  In  addition  the  compressive  load  causes  a  compressive  stress 

<7  =  25000  psi 
x 

i.e.  cr  total  =  26262  psi. 
x 
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0.4  0.5  0.75  1.0  1.5  2.0  2.5 

Aspect  Ratio  a/b 

Figure  III-l.  Maximum  Deflection  (W  )  °f  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (  A  T)  Between  the  Faces  and  Subjected  to  a  Compression  77  =  0.0.  All 

Edges  Simply  Supported 
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Figure  IH-4.  Maximum  Deflection  of  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (AT)  Between  the  Faces  and  Subjected  to  a  Compression  rj  =  -0.15 

All  Edges  Simply  Supported 
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Figure  III-5.  Maximum  Deflection  (W  )  of  a  Rectangular  Sandwich  Panel  Due  to  a  Uniform  Normal 
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Figure  III- 6.  Maximum  Deflection  (Wm)  of  a  Rectangular  Sandwich  Panel  Due  to  a  Uniform  Normal 
Pressure  (p>  and  Subjected  to  a  Compression  -rjx  =  -0.05.  All  Edges  Simply  Supported 
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Figure  III-7.  Maximum  Deflection  (Wm)  of  a  Rectangular  Sandwich  Panel  Due  to  a  Uniform  Normal 
Pressure  ( p>  and  Subjected  to  a  Compression  tj  -  -0.10.  All  Edges  Simply  Supported 
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Figure  III-8.  Maximum  Deflection  (Wm)  of  a  Rectangular  Sandwich  Panel  Due  to  a  Uniform  Normal 
Pressure  (p>  and  Subjected  to  a  Compression  77  =  -0.15.  All  Edges  Simply  Supported 
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Figure  III—  9.  Maximum  Bending  Moment  In  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (  A  T)  Between  the  Faces  and  Subjected  to  a  Compression  77 x  =  0.0. 

All  Edges  Simply  Supported 
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Figure  III-10.  Maximum  Bending  Moment  In  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (  AT)  Between  the  Faces  and  Subjected  to  a  Compression  77  =  -0.05 

All  Edges  Simply  Supported 
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Figure  in- 11.  Maximum  Bending  Moment  in  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature  Difference 
(AT)  Between  the  Faces  and  Subjected  to  a  Compression  17^  =  -0.10.  All  Edges  Simply  Supported 
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Figure  III  12.  Maximum  Bending  Moment  in  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature  Difference 
(A  T)  Between  the  Faces  and  Subjected  to  a  Compression  tj  =  -0.15.  All  Edges  Simply  Supported 
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Figure  III- 15.  Maximum  Bending  Moment  in  a  Rectangular 
(p)  and  Subjected  to  a  Compression  rj  = 


Figure  III- 16.  Maximum  Bending  Moment  in  a  Rectangular  Sandwich  Panel  Due  to  a  Uniform  Normal  Pressure 
(p)  and  Subjected  to  a  Compression  ^  =  -0.15.  All  Edges  Simply  Supported 


-0.0 


Jj 


Ml 


jWi; 


Midplane  Loading  t)  x  =  Nx/Dq 

-18.  Maximum  Defleetion  of  a  Rectangular  Sandwieh  Panel  Due  to  a 
Uniform  Normal  Pressure(p)  and  Subjected  to  Midplane 
Compression  All  Edges  Simply  Supported. 
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Figure  III-19.  Maximum  Deflection  (W  )  of  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (AT)  Between  the  faces.  Two  Edges  Simply  Supported  and  Two 

Edges  Built-In. 
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Figure  III-20.  Maximum  Deflection  (Wm)  of  a  Rectangular  Sandwich  Panel  Due  to  Uniform 
Normal  Pressure(p).  Two  Edges  Simply  Supported  and  Two  Edges  Built-In 
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Figure  III- 21.  Bending  Moments  at  Center  of  a  Rectangular  Sandwich  Panel  Due  to  a  Temperature 
Difference  (AT)  Between  the  Faces.  Two  Edges  Simply  Supported  and  Two  Edges  Built-In 


o  Edges  Simply  Supported  and  Two  Edges  Built- 
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CHAPTER  IV 

SANDWICH  CYLINDER  INSTABILITY  UNDER  CIRCUMFERENTIALLY 
VARYING  AXIAL  STRESS 


This  chapter  is  concerned  with  the  prediction  of  the  elastic  instability  of  an 
axially  loaded  sandwich  cylinder  of  mean  radius  r  and  core  thickness  h,  with 
isotropic  faces  of  equal  thickness  t  (Figure  IV-1).  The  faces  possess  only  inplane 
stiffness  and  have  a  modulus  of  elasticity  E,  while  the  core  possesses  only  a 
transverse  shear  stiffness  with  a  modulus  of  rigidity  Gc>  The  axial  load  per  unit 
circumferential  length,  Nx,  varies  in  the  circumferential  (y)  direction  and  at  any 
point  the  stress  in  each  of  the  two  faces  is  assumed  equal  at  a  value  av  =  jjx 

X  2t 

Conditions  of  circumferentially  varying  stress,  as  shown  here,  may  result 
from  the  combination  of  an  applied  axial  load  and  bending  moment,  or  as  a  result 
of  temperatures  which  vary  around  the  circumference  at  the  cylinder.  Trans¬ 
formation  of  the  applied  loads  or  temperature  profile  into  the  required  stress 
distribution  is  assumed  to  have  been  accomplished  by  the  user  of  the  present  data 
by  means  of  the  appropriate  stress  analysis  technique. 


The  problem  of  the  elastic  instability  of  an  isotropic  thin-walled  cylinder 
has  been  solved  in  References  9  and  10  on  the  basis  of  small  deflection  theory. 
Both  references  concluded  that  instability  is  reached  when  the  maximum  axial 
stress  has  the  value 


a  ~  0.6  Et 

cr  - 

r 


(IV-1) 


regardless  of  the  nature  of  the  circumferential  variation  of  <xx.  (The  compression 
zone  within  which  the  maximum  value  crcr  occurs  must  at  least  extend  over  the 
wavelength  of  a  circumferential  buckle,  however).  "t  is  the  wall  thickness  of  the 
isotropic  shell.  The  result  of  Equation  IV-1  is, of  course,  the  small  deflection 
theory  solution  for  uniform  compression.  In  view  of  these  results  it  would  appear 
reasonable  to  expect  that  for  the  case  of  a  sandwich  cylinder  the  same  conclusions 
would  prevail,  i.e.,  the  instability  stress  is  given  by  the  solution  for  uniform  axial 
compression,  regardless  of  how  nonuniformly  the  stress  might  vary. 


Small  deflection  theory  solutions  for  the  instability  at  sandwich  cylinders 
under  uniform  axial  compression,  based  on  infinite  series  techniques,  have  been 
presented  in  references  11,  12  and  13.  These  results  are  useful  for  comparison 
purposes  but  to  develop  solutions  for  nonuniform  stress  it  has  been  found  more 
convenient  to  extend  the  finite  difference  technique  used  in  Reference  9  to  the 
present  case.  Details  of  this  extension  are  given  in  Appendix  C. 
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Based  on  the  theory  described  in  Appendix  C  a  computer  program  for  the 
prediction  of  critical  stress  under  the  subject  conditions  was  coded  and  employed 
in  the  development  of  parametric  results.  For  sandwich  cylinders,  two  stiffness 
parameters  must  be  considered: 


D 


c 


(IV- 2) 


/  Et  \ 

and(- — -j.  The  chosen  values  of  Dc  range  from  66.7  to  2500,  while 
(jcr 


Et 

Gcr 


ranges  from 


0  to  approximately  0.95.  The  upper  limits  of  both  parameters  are  governed  by  a 
wrinkling  failure  characterized  by  a  shear  instability  of  the  core.  As  in  Chapter  II, 
this  instability  is  governed  by  the  condition 


2  ^cr1 
(h+t)Gc 


=  1 


(IV- 3) 


The  specific  load  conditions  studied  are  those  which  produce  a  linear  variation 
of  stress,  described  by  the  ratio  (S)  of  the  crown  stress  (  c tXq)  to  the  stress  at  the 
bottom  of  the  cylinder  (  cr^),  Ee., 


S  = 


(IV-4) 


Results  were  obtained  for  S  =  1.0,  0.5,  0,  -0.5  and  -1.0  for  various  combinations,  of 
the  two  governing  parameters. 


Selected  results  are  presented  in  Table  IV-1.  As  anticipated  the  critical 
stresses  for  nonuniform  stress  states  are  effectively  equal  to  the  critical  stress 
for  uniform  axial  compression  (S  =  1.0).  The  negligible  differences  between  the 
results  for  a  given  stiffness  condition  may  be  the  result  of  the  differences  in  stress 
distribution,  but  it  is  also  possible  that  the  discrepancies  are  due  entirely  to  numerical 
error.  This  question  cannot  be  resolved  by  use  of  the  method  adopted  for  solution  of 
the  problem. 


Results  are  plotted  in  Figures  IV-1  through  IV-3.  Figure  IV-1  is  a  carpet 
plot  wherein  both  stiffness  parameters  play  a  role  in  the  definition  of  the  critical 
stress.  If  a  conventional  plotting  procedure  and  slightly  altered  scales  are  adoped, 
as  in  Figure  IV-2,  the  parameter  Dc  loses  its  significance  and  all  results  can  be 
approximated  by  a  single  line.  This  manner  of  representation  was  also  adopted  in 
References  11,  12,  and  13.  The  curve  presented  in  Figure  IV-2,  based  on  an 
Equation  from  Reference  13,  is  given  by 
E(h  t) 


•H 


(IV— 5) 
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TABLE  IV-1 

BUCKLING  COEFFICIENTS  -  SANDWICH  CYLINDER 
UNDER  NONUNIFORM  AXIAL  COMPRESSION 


Stiffness 

Parameters 

2t  CT 
cr 

Critical  Stress  Parameter  -  Q  (computed) 

q 

D 

c 

Et/G  r 
c 

m 

m 

S  =  0 

S  =  0.5 

S  =  1.0* 

1000 

13.49 

0.2672 

Q.2669 

0.2665 

0.2660 

0.2630 

200 

30.17 

0.5536 

0.5511 

0.5476 

0.5433 

0.5345 

66.7 

45.24 

0.7699 

0.7630 

0.7555 

0.7450 

0.7243 

*  Uniform  axial  compression 
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The  present  results,  as  shown  by  individual  points  on  Figure  IV-2  are  in  very  good 
agreement  with  the  above  curve. 

For  the  isotropie  cylinder,  results  obtained  from  small  defleetion  theory 
analyses  are  grossly  in  error  except  when  a  sufficiently  high  internal  pressurization 
is  applied,  the  attainable  critical  stress  being  extremely  sensitive  to  initial 
imperfections  (e.g.,  out  of  roundness,  etc).  It  is  therefore  to  be  expected  that  sand¬ 
wich  cylinder  critical  stresses  will  be  in  closer  agreement  with  small  deflection 
theory  predictions  since  the  initial  imperfection  effeets,  whieh  are  a  function  of  the 
radial  imperfection  magnitude-to-total  wall  thickness  ratio,  will  be  very  small  in 
carefully  fabricated  sandwich  cylinder. 

Experimental  results  presented  by  Cunningham  and  Jaeobson  (Reference  24) 
lend  eredenee  to  the  above  hypothesis,  On  the  other  hand  the  test  data  of  Norris  and 
Kuenzi  (Reference  25)  support  the  view  that  largc-deflection  formulations  must  be 
used  for  the  development  of  design  data.  A  design  curve  from  Reference  26,  based 
on  large  deflection  theory,  is  reproduced  in  Figure  IV-2.  Note  the  large  differences 
between  the  predictions  based  on  large  and  small  defleetion  theory.  Unquestionably, 
a  rational  basis  for  the  design  of  sandwich  cylinders  requires  a  correlation  and 
critical  evaluation  of  existing  test  data,  with  the  performance  of  additional  tests  in 


work  should  be  implemented  by  a  theoretical  study  whieh  treats  instability  under 
nonuniform  stress  on  the  basis  of  large  deflection  theory.  A  study  of  this  type  is 
beyond  the  scope  of  the  present  effort. 

As  a  final  representation,  the  variation  of  the  wavelength  of  the  buekle  in  the 
axial  direction  is  plotted  parametrically  in  Figure  IV-3.  As  in  the  case  of  the 
eigen  values,  the  computed  wavelengths  were  found  to  be  essentially  independent  of 
the  shape  of  the  stress  distribution.  The  wavelength  to  radius  ratios  are  relatively 
small;  thus,  the  present  results  should  apply  to  short  as  well  as  to  long  cylinders. 

It  should  be  noted  that  these  eomputed  wavelengths  are  in  elose  agreement  with 
results  presented  elsewhere  (Reference  12). 

Illustrative  examples 

(1)  "Rigid  eore" 

Determine  the  maximum  stress  for  buekling  (  <j  )  for  an  axially  loaded 
sandwieh  cylinder  possessing  the  following  properties: 

Radius  (r)  :  18  in. 

Core  Depth  (h)  :  0.125  in. 

Faec  Thiekncss  (t):  0.010  in. 

Moduli  E  =  10.5  x  10^  psi 
Ge  =  40,000  psi 
fJL  =  1/3 
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Thus,  D 


2(l-yu. 2)  G^r2  =  1660 
Et  (h+t) 


— —  =  0.146 
G  r 
c 

From  Figure  IV- 1:  k=0,96 

,  „  (h+t)  0.96  x  10.5  x  106  x  0.135 

Hence  <J  =  k  E  - — -  =  - — - 

cr  r  18 


75,500  psi 


If  the  material  in  question  is  7075-T6  aluminum,  as  was  the  case  for 
uniformly  compressed  cylinders  of  the  same  proportions  tested  by  Eakin  (Reference 
27)  and  discussed  by  Cunningham  and  Jacobson,  this  stress  is  well  into  the  inelastic 
range.  By  using  the  dotted  curve  of  Figure  IV-2  and  the  appropriate  plasticity 
reduction  factor,  Cunningham  and  Jacobson  found  a  small-deflection  theory  critical 
stress  of  57,000  psi.  The  test  specimen  failed  at  61,000  psi. 


(2)  "Soft  Core" 

Assume  all  of  the  above  properties  and  proportions  remain  the  same 
except  that  the  core  now  has  a  shear  modulus  Gc  =  1450  psi  and  h  =  0.1420  in. 

The  parameter  Et/Gcr  will  then  be  equal  to  the  value  prevailing  for  tests  described 
in  Reference  12.  Thus, 

Et  10.5  x  106  x  0.0100  ,  _ 

- =  -  =  4.02 

Gcr  1450  x  18 

From  Figure  IV-2,  it  is  seen  that  failure  is  governed  by  wrinkling  associated  with 
shear  instability  of  the  core.  Hence 
(h+t)Gc 

^cr  =~  =  10,300  psi 


This  compares  with  a  stress,  computed  from  experimental  data, of  11,300  psi. 
(Specimen  1424A,  Reference  12). 
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Figure  IV- 1.  Variation  of  Critical  Reference  Stress  with  Two  Stiffness  Parameters 
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CHAPTER  V 

THERMAL  STRESSES  AND  BENDING  MOMENTS  IN  A  HEATED  SANDWICH 
CYLINDER  SUPPORTED  BY  RIGID  BULKHEADS 


This  ehapter  presents  data  for  the  determination  of  the  thermal  stresses  and 
displacements  in  a  sandwich  cylinder  of  length  2e  supported  at  eaeh  end  and  sub¬ 
jected  to  a  temperature  ehange  (  A  T)from  the  stress-free  state.  The. cross-sec¬ 
tional  properties  and  the  material  property  characteristics  are  identical  to  those  of 
the  sandwieh  ey Under  of  the  preeeding  chapter  (see  Figure  V-l).  Both  fixed  and 
simple  support  conditions  are  treated.  The  fixed  support  condition  applies  to  a 
cylinder  continuous  over  many  supports,  where,  due  to  symmetry,  each  bay  is 
effectively  built-in  at  the  bulkheads. 

A  detailed  development  of  the  pertinent  formulations  is  presented  in  Appendix 
D.  There,  it  is  shown  that  the  problem  is  similar  to  that  of  a  beam  on  an  elastie 
foundation.  The  radial  displacements  (w)  for  the  simple  support  condition  are  found 
to  be 


rq  A  T 
(m22-m12) 


jm22-2) 


2  eosh  miHY  ,  2  „  eosh  moH*  ,  2 
-2)  - L2l  -  (m  -2)  - (m  -  m] 

eosh  m^Hj  eosh  n^H 


where  m. 


1  +  \/l  - 


Hi 


EtDs  Et 

'’2D  2  D  D 
q  ucuc\ 


‘1  r  v  D 


e  /Et 


The  hoop  stress  (  cr^  )  midway  between  the  supports  (at  x  =  0),  which  is  a  ring 
compression  or  tension  (i.e.,  it  nas  an  equal  value  in  both  the  inner  and  outer  faees) 
is  given  by 


E  a  AT 


m-,2-2  m„2-2 


e  (m^-nij2)  Lcosh  m  Hj  eosh  m^II^ 


while  the  longitudinal  bending  moment  (Mx  )  at  x  =  0  is  ealeulable  from 
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My  H^D  raAT  (mi2-2)  (m22-2) 

Ac  =  a  q  - - - 

(m22_mi2) 

Formulas  for  the  center  deflection,  the  hoop  stress  at  the  center,  and  for  the 
longitudinal  bending  moments  at  the  center  and  supports  of  the  cylinder  are  given  by 
Equations  D-9  through  D-12  of  Appendix  D. 

The  hoop  stresses  and  bending  moments  for  both  support  conditions  were  com¬ 
puted  for  realistic  ranges  of  the  pertinent  parameters  (Hj  and  Hj)  and  are  plotted  in 
the  form  of  nondimensionalized  design  charts  in  Figures  V-l  through  V-5.  These 
figures  show  the  characteristic  reversal  of  stresses  associated  with  the  beam  on  an 
elastic  foundation  type  of  problem. 

Illustrative  Examples 

(1)  "Isotropic"  cylinder 

As  a  check  on  the  accuracy  of  the  plotted  data  the  case  of  an  isotropic 
cylinder  with  simply  supported  ends,  examined  by  Hoff  (Reference  3)  will  be 
treated  (there  are  no  alternative  solutions  published  for  the  sandwich  problem). 

In  Hoff's  paper,  r  =  10",  C  =  1.57",  T  =  0.0331"  and  E  =  29  x  103  psi.  To  achieve 

a  comparison  with  the  present  paper  an  equivalent  sandwich  must  be  defined.  It 
has  been  suggested  that  a  suitable  criterion  is  that  the  radii  of  gyration  of  the 
sandwich  and  the  isotropic  plate  shall  be  equal.  Also  the  total  cross  sectional 
areas  must  be  equal. 


cosh  m2Hj  cosh  m^Hj 


(V-8) 


Thus,  if  t  is  the  thickness  of  isotropic  sheet,  t  is  the  thickness  of  each  sand¬ 
wich  face  and  h'  is  the  effective  depth,  (i.e.  distance  between  centroid  of  the  faces). 

th'2 


—  f 
then  2t  =  t  and  = 


2.2t 


i.e.  t=-~  and  h' =  7= 
2  V  3 


Also  a  core  shear  stiffness  Dqinust  be  defined.  However,  in  the  case  of  an 
isotropic  plate  the  shear  deformation  has  negligible  influence  and  a  large  value 
for  Gc  may  be  chosen  arbitrarily.  Using  the  formula  Dq  =  h'  Gc  and  with  an 
arbitrary  choice  of  Gc  =  1.1  x  103  psi  one  obtains  Dq  =  2.1  x  104.  The  bending 
stiffness  is  computed  as 


Ds 


96.307 


hence 


EtD 


H 


d  rV 
q 


and  H^ 


0.742 
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From  Figure  V-2  at  =  10_,i  and  =  0.742  the  hoop  stress  at  center 
ar^  =  0.057  E  a  AT 

In  Reference  3  the  stress  at  the  center  of  the  cylinder  is  obtained  from  the 
expression 


cosh  /3  cos  /3 

$  sinh2  /3  +  cos2/3  J 

I  ac2  I/O 

where  /3  -  0.643  =  3.51 


E  a  At 


The  stress,  as  dervied  from  the  above  formula,  is  then 
o-i  =  0.056  Ea  AT 

Thus,  there  is  excellent  agreement  between  the  graphical  and  analytical  results. 

(2)  Sandwich  Cylinder 

A  long  sandwich  cylinder  of  70  inches  mean  radius  is  ring  stiffened  at 
intervals  of  48  inches.  The  core  thickness  is  1.0  inch,  the  face  thicknesses  0.030 
inch.  Both  the  core  and  the  faces  are  composed  of  17-7PH  stainless  steel,  E  =  27  x 
10®  psi,  d  =  6.1  x  10"®  /°F,  Gc  =  100,000  psi.  The  cylinder  is  uniformly  heated 
to  a  temperature  600°F  in  excess  of  the  ring  temperature.  Determine  the  resulting 
hoop  stresses  and  longitudinal  bending  moments. 


Assume  the  ring  stiffeners  arc  infinitely  rigid  and  that  in  any  particular 
bay  the  cylinder  can  be  assumed  to  be  fixed-supported  at  its  ends. 


D  =  J  ^  =  1070 

c  Et(h+t) 


D  =  (h+t)  G  =  103,000 
q  c 


II  =  =  7.35  x  10  ^ 

d  DcDq 


Vf/ij  =°-961 


From  Figure  V-3  by  interpolation  for  and  the  hoop  stress  is  given  by 


o-^j  =  0.060  E  a  AT  =  5930  psi 


From  Figures  V-4  and  V-5  the  bending  moments  at  the  center  and  ends  arc  found  to 
be 


At  center  M  =  0.0152  D  ra  A  T  =  401  lb  in./in. 
xc  q 

At  ends  Mv  =0.12D  ra  AT  =  3165  lb  in./in. 
xe  q 
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Figure  V  1.  Central  Hoop  Stress  in  Sandwich  Cylinder  Simply  Supported  at  Each  End  with 
Rigid  Bulkheads  Subjected  to  Temperature  Increase  AT 
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CHAPTER  VI 

INSTABILITY  OF  HEATED  SANDWICH  CYLINDERS  SUPPORTED 
BY  RIGID  BULKHEADS 


The  data  presented  in  this  Chapter  permit  the  determination  of  the  temper¬ 
ature  change  which  will  produce  the  elastic  instability  of  sandwich  cylinders 
supported  by  rigid  bulkheads.  The  conditions  of  analysis  are  similar  to  those  of 
the  preceding  Chapter,  where  the  stresses  introduced  by  the  temperature  change 
were  determined.  In  determining  the  stresses  only  simply  supported  and  fully 
fixed  cylinders  were  considered.  However,  in  considering  the  instability  of  the 
cylinders  under  hoop  stress  alone,  it  is  necessary  to  define  three  conditions  of 
analysis. 

a)  Fixed  support  condition.  In  this  case  the  cylinder  consists  of  a  single 
bay  fully  fixed  against  rotation  by  a  rigid  bulkhead  at  each  end. 

b)  Simple  support  condition — single  bay.  Again  a  single  bay  cylinder  is 
considered  but  the  rigid  bulkheads  no  longer  restrain  the  rotation  of  the  ends. 

c)  Simple  support  condition--many  bays.  The  cylinder  is  taken  to  be 
continuous  over  many  bays,  as  in  a  fuselage.  The  bulkheads  restrain  radial  dis¬ 
placement  but  provide  no  restriction  against  rotation. 

In  the  case  of  a  long  multi-bay  cylinder  (Case  c))  the  thermal  stress  dis¬ 
tribution  in  each  bay  corresponds  to  the  fully  fixed  case.  The  fixing  moment  is 
produced  by  symmetry  at  the  supports.  Although  one  possible  buckling  mode  could 
be  the  fully  fixed  one,  as  discussed  in  Section  2  of  Appendix  E,  where  each  bay  is 
fully  fixed  due  to  symmetry  at  the  bulkheads,  buckling  could  also  occur  in  the 

simply  supported  mode  as  shown  in  the 
adjacent  sketch.  In  this  latter  case 
each  bay  behaves  as  if  it  were  simply 
supported  at  the  bulkheads  and  since 
this  mode  will  involve  less  bending  of 
the  faces  than  the  fully  fixed  mode 
(see  sketch),  the  simply  supported 
mode  will  predominate. 
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Thus,  for  a  multi-bay  cylinder  of  this  type  the  stress  distribution  must  be 
determined  for  the  fully  fixed  case  and  this  stress  distribution  used  in  the  simply 
supported  stability  analysis  to  determine  the  critical  temperature  Tcr.  This 
method  has  been  used  to  develop  curves  for  the  determination  of  the  critical 
temperature  of  a  simply  supported  multi-bay  sandwich  cylinder.  These  curves  are 
presented  in  Figure  VI-2. 

Figure  VI-1  provides  the  means  for  determining  the  critical  temperature  for 
a  simply  supported  single  bay  (case  b)). 


The  pertinent  parameters  are  in  the  same  range  of  values  as  were  employed 
for  the  design  curves  of  Chapter  V. 

t  E  a  AT 

cr 

Both  sets  of  curves  have  been  terminated  at  an  upper  value - — -  =  O.E 

Dq 

corresponding  to  the  shear  instability  failure  of  the  core  noted  previously.  This 
instability  was  characterized  by  the  relationship 


At  the  restrained  ends  of  the  cylinder  the  hoop  stress  (( )  equals  E<X  A  T  . 
Thus  01 

2t  E  <X  A  T 

cr 

D  =1'° 

q 


or 

t  Ea  AT 

q 

The  significance  of  the  present  data  with  respect  to  published  isotropie 
cylinder  solutions  is  discussed  at  length  in  the  first  of  the  two  following  illus¬ 
trative  examples. 

(1)  Isotropic  Cylinder 

6  6 

A  steel  (E  =  29  x  10  psi,  <X  =6.1  x  10  /°F)  isotropic  cylinder  of 
radius  15  in.,  wall  thickness  0.022  in.,  is  simply  supported  at  rigid  bulkheads  at 
intervals  of  5  in.  The  eritieal  temperature  is  to  be  determined  for  each  of  two 
cases: 
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For  this  cylinder 


2 

(2c) 

r  T 


25 

15  x  0.022 


=  75 


In  Figure  1  of  Reference  18,  Anderson  presents  the  critical  temperature  for  an 
(2c)2 

isotropic  cylinder  ^  =  75.  The  corresponding  values  of  (OfAT^) 

as  taken  from  the  curve,  are  5.24  and  2.08  respectively.  The  actual  critical 
temperatures  are 


a) 


cr 


0.121  D 
_ c 

t  EC?' 


=  1245  F 


b) 


cr 


0.049  D 
_ c 

t  E0C 


=  504  F 


(2)  Sandwich  Cylinder 


A  sandwich  cylinder  of  95  in.  mean  radius  is  ring  stiffened  at  intervals 
48  in.  The  core  thickness  is  1.0  in.,  the  face  thicknesses  0.030  in.  Both  core  and 

faces  are  composed  of  17-7  PH  stainless  steel,  E  =  30  x  106  psi,  Ot  =  6.1  x  10  ^ 
in./in.°F,  G  =  70,000  psi.  The  critical  temperatures  are  again  to  be  determined 
for  the  simply  supported  single  and  multi-bay  cases. 


Then  D  =  h  G  =  7.21  x  104 
q  c 


E  t  E> 

and  Hd  =  -j-f  =  10- 

r  D 

q 


and  D  = 
s 


and  II.  - 


E  th 

2(1-/ 


=  5.246  x  10 


Etc 


2‘ 


r2  D 


1/2 


=  0.893 


q  J 


From  Figure  VI-1  for  =  10“^  and  Hj  =  0.893  the  point  lies  above  the 
shear  instability  boundary.  Thus  the  critical  temperature  for  the  single  bay 
cylinder  is  given  by 


and 


t  E0C  AT 

cr  „  _ 

5  =0-5 

q 

0.5  D 

AT  -  — n~.  6566°  F 

cr  t  E  0C 
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From  Figure  VI-2  for  the  relevant  and  H[  for  a  multi-bay  cylinder 


t  E  Of  AT 

cr 

- -  .  0.4 

q 

and  AT  =  5253°F 

cr 

It  can  be  seen  from  the  above  examples  that  this  type  of  instability  is 
unlikely  to  occur  in  sandwich  cylinders  of  normal  proportions.  It  is,  however,  a 
possible  mode  of  failure  in  cylinders  of  very  large  radii  and  low  wall  thickness. 
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CHAPTER  Vn 

SANDWICH  CYLINDER  STRESSES  DUE  TO  INTERNAL  AND  EXTERNAL 
PRESSURE  AND  RADIAL  TEMPERATURE  GRADIENTS 


Design  curves  presented  in  this  section  permit  the  determination  of  the 
stresses  in  long  sandwich  cylinders  subjected  to  radial  temperature  gradients  and 
to  internal  and  external  pressure.  As  in  the  preceding  chapter,  the  cylinder  consists 
of  two  concentric  thin  cylinders  of  identical  materials,  with  radii  rQ  and  r.  and 
thicknesses  t,  separated  by  a  honeycomb  type  core  of  thickness  h  =  rQ-r-  (see 
Figure  VII-1).  The  cylinder  is  heated  such  that  each  skin  has  a  constant  tempera¬ 
ture  change  from  the  stress-free  state  (T0  and  T^,  respectively)  with  a  linear  vari¬ 
ation  across  the  core,  while  the  internal  and  external  pressures  are  p0  and  p^. 

The  theoretical  basis  of  the  present  chapter  is  detailed  in  Appendix  F,  where 
the  desired  stress  equations  are  derived  by  first  considering  the  separate  deforma¬ 
tions  of  the  three  components  -  -  the  two  skins  and  the  core — under  the  actions  of 
external  and  interface  loading  and  by  then  applying  the  conditions  of  interface  com¬ 
patibility.  The  isotropic  faces  are  analyzed  using  plane  stress  relationships.  The 
core  is  assumed  to  possess  a  modulus  of  elasticity  in  the  radial  direction  (Er)  and 
a  coefficient  of  thermal  expansion  in  the  radial  direction  equal  to  that  of  the  skins. 
Consistent  with  the  common  properties  of  honeycomb  construction,  the  core  has 
neither  stiffness  nor  a  coefficient  of  thermal  expansion  in  the  axial  or  circumferen¬ 
tial  directions.  As  shown  in  the  illustrative  example,  however,  homogeneous 
(rather  than  honeycomb)  cores,  of  the  type  often  employed  in  sandwich  construction, 
can  be  analyzed  with  sufficient  accuracy  by  use  of  the  present  curves. 


It  should  be  noted  that  the  ends  of  the  cylinder  are  assumed  to  be  free  to  move 
in  the  axial  direction  and  axial  stresses  due  to  end  closure  of  the  cylinder  are  not 
included.  The  latter,  which  are  approximately  equal  tol-L-i,  can  be  superimposed 
on  the  stresses  obtained  from  the  present  curves.  "End  effect"  stresses,  if  calcu¬ 
lable,  can  also  be  superimposed. 


Due  to  the  imposed  load  and  temperature  conditions,  axial  and  longitudinal 
stresses  are  produced  in  the  faces  and  radial  direct  stresses  are  introduced  into  the 
core.  There  are  no  direct  axial  or  circumferential  stresses  in  the  core  due  to  the 
absence  of  core  stiffness  in  these  directions.  Design  curves  for  the  axial  and  cir¬ 
cumferential  stresses  at  the  outside  surface  due  to  external  pressure,  internal 
pressure,  and  temperature  are  presented  in  Figures  VII-1  through  VII-6.  The  rela¬ 
ted  formulas  for  stress  are  shown  in  each  figure;  these  formulas  show  that  the 
dimensionless  stress  parameters  (e.g.  ag  T/Ed  (Tq-T^))  depend  only  upon 

o 

R  =  r0/rj  and  Et/E^.  For  shells  of  current  interest  and  importance,  appropri¬ 
ate  values  of  R  lie  between  1.0025  and  1.10  while  the  ratio  Et/E^  takes  on  values 
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in  the  range  of  0.01  to  50.0.  These  ranges  were  employed  in  the  development  of  the 
subject  design  curves. 

Values  for  the  inner  surface  and  radial  stresses  are  easily  obtained  from  the 
corresponding  outer  surface  stresses  through  simple  algebraic  operations  and  are 
therefore  not  plotted.  The  pertinent  relationships  are:  (Equations  VII- 1  to  VII-6 
appear  in  Figures  VII- 1  to  VII-6). 


a.  Stresses  Due  to  Radial  Temperature  Gradient 


T 

T 

9  i 

=  ~a9  o 

T 

T 

=  -  R  a 

x , 

i 

xo 

(VII-7) 
(VII- 8) 

(VII-9) 


b.  Stresses  Due  to  Internal  Pressure 


(VII- 10) 

(VII-11) 

(VII- 12) 


c.  Stresses  Due  to  External  Pressure 


(VII— 13) 

(V 11—14) 

(VII- 15) 


Illustrative  Example 

Given  a  sandwich  cylinder  with  l'i  =  100.01  in.,  r0=  101.03  in t  =  0.020  in., 
subjected  to  a  radial  temperature  gradient  TQ  -  T.  =  100  °F,  and  pressures  p^  =  10 

psi  and  P0  -  5  psi.  The  skin  material  is  aluminum  alloy  with  E  =  10.5  x  10^  psi  and 
/l  =  0.3  ,  while  the  core  is  aluminum  honeycomb  with  Er  =  1.5  x  10^  psi.  For 
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Radial  stress  in  core 

T  Pi  Po 

a  =  +1.8  -  4.9  +  2.5  =  -0.5  psi 
r 

Similarly,  the  inner  skin  stresses  are  found  to  be 

a  =  -21,900  psi  and  &  =  +9010  psi. 

9  i  Xi 

(19) 

Yao  lias  analyzed  the  case  of  a  sandwich  cylinder  under  a  radial 
thermal  gradient  alone,  where  the  core  possesses  the  properties  of  a  homogeneous 
material.  Otherwise,  the  problem  data  are  identical  to  those  given  above.  The 
stresses  found  by  Yao  are  as  follows: 

T 

cr  =  -0.730  E  a  (T0-Tj)  =  -9200  psi 

V  O 

T 

a  9  .  -  0.085  E  a  (Tq-T.)  =  8630  psi 
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cr  J  =  1.37  x  10"4  E  a  (T^Tj)  =  1.7  psi 

These  values  compare  closely  with  the  present  results  of  -9000,  9000,  and  1.8  psi, 
respectively. 
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Figure  VII-2.  Hoop  Thermal  Stress  in 
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Sandwich  Cylinder  Due  to  Internal  Pressure  ( 
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Figure  VII-4.  Hoop  Stress  in  Outer  Face  of  Sandwich  Cylinder  Due  to  Internal  Pressun 
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APPENDIX  A 

BUCKLING  OF  SANDWICH  PANELS  UNDER 
NONUNIFORM  STRESS 


The  governing  differential  equation  for  the  stability  of  a  sandwich  plate 
composed  of  equal  thickness  isotropic  faces  possessing  only  stiffness  in  direct 
stress  and  an  isotropic  core  which  possesses  only  shear  stiffness  (see  Figure  II-l 
in  text)  can  be  written  in  the  form  (Reference  20) 

Ds  V4W  +  (1  3-  V  2)  2t  cr  =  0  (A-l) 

q  x  j 


where  Dg  is  the  panel  bending  stiffness 
D„  is  the  core  shear  stiffness 


4  2  2 

V  =  V  V 

To  solve  Equation  (A-l),  the  displacement  is  first  assumed  as 
w  =  W  sin  X  x 


where  W  is  a  function  of  y  only  and  X  is  a  parameter  involving  the  buckle 
wavelength  in  the  longitudinal  (x)  direction.  Using  Equation  (A-2)  in  (A-l)  and 
introducing  a  nondimensionalization  of  the  stress  distribution  through  use  of 


cr  x 

P  y  - - ,  one  obtains 

a  cr 


X  W  -2  X  W  +  W 


'1  2  ta~c 


/  °s  *  \  °s  n 
i  +  - — )p  w -—(p  W)  =c 
\  /  y  d  y  J 


where  the  primes  on  W  and  (  p  ^W)  indicate  derivatives  with  respect  to  y.  Note 
that  the  introduction  of  p  serves  to  establish  a  single  parameter  ( crcr)  which 
characterizes  instability  in  the  presence  of  nonuniformly  distributed  stresses. 

Finite  differences  are  used  to  reduce  this  equation  to  algebraic  form.  With 
9  intervals  (see  Figure  II- 1)  and  second  order  finite  differences  one  obtains  the 
following  evaluation  of  Equation  (A-3)  at  an  arbitrary  point  k 
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-2  Wk-3  +  Wk+3 


+  W,  „+W.  n 
k-2  k+2 


(24  +-  A  ) 


-  W  +  W  I  (78  +|f  A  )  +  W.  (12  +|2.  A  +  12  (— )  2)1 
k-1  k+1  J  '  81  ’  k  '  27  \8l/  l\ 

_  IT  A  /a  w  -16A  w  *p  w  (30+i-A  +  i-DP  ) 
81  k-2  k-2  k-1  k-1  k  k  27  27 


16  p  W  p  w 
k+1  k+1  +  k+2  k+2 

2t  a  D  b2 

where  FI  =  cr  ,  Dp=  — —  and  A  =  (  X  b) 

q  s 


The  boundary  conditions  on  W  are: 

For  simple  support 

WQ  =  W9  =  0,  W_:  =  -W:  ,  W_2  =  -W2  ,  W1Q  =  -Wg  and  Wn  =  -  W? 

For  fixed  support 

WQ  =  W9  =  0,  W_j  -  Wp  W_2  =  W2  ,  W1Q  =  Wg  and  Wu  =  W? 

The  boundary  conditions  on  p  must  be  based  on  an  extrapolation  of  the  stress  dis¬ 
tribution  beyond  the  edges  of  the  plate.  (Although  there  is  no  loading  outside  of  the 
edges  of  the  plate  (i.e. ,  in  actuality  p  j  =  0  for  i  <0  and  i  >  9)  the  use  of  finite 
differences  requires  the  introduction  of  fictitious  exterior  loadings  if  the  stress 
variation  in  the  intervals  adjacent  to  the  edges  are  to  be  represented  properly). 

With  use  of  the  boundary  conditions,  Equation  (A-4)  can  be  adapted  to  each  of 
the  eight  internal  points,  obtaining  thereby  eight  simultaneous  equations  with 
coefficients  W^,  — Wg.  These  equations  can  be  written  in'  matrix  form  as  follows 


[a]  -  n a[b] 


W  =0 


The  matrices  J^'aJ  and  bJ  are  presented  in  Figures  A-l  and  A-2,  respectively. 
Both  simply  supported  and  fully  fixed  edges  are  included. 

The  condition  for  buckling  is  that  the  determinant  of  Equation  (A-5)  shall 
vanish.  A  more  attractive  solution  procedure  is  through  matrix  iteration,  however, 
since  direct  evaluation  of  the  determinant  would  lead  to  a  somewhat  complex 
polynomial  of  eighth  order  in  IT  .  To  permit  matrix  iteration,  Equation  (A-5)  is 
rearranged  as  follows: 

1  ,  . 


M  =  [a]  -1  [b]  {w} 
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Then,  iteration  of  Equation  (A-6)  is  performed,  starting  with 
vector  {w  }  ,  until  convergence  on  the  eigenvector  |w|  cr 


is  achieved. 


an  initially  assumed 
and  eigenvalue 


For  this  type  of  problem  it  has  been  found,  under  certain  conditions,  that  an 
oscillation  will  occur  in  the  iterative  process.  The  process  converges  on  two  dif¬ 
ferent  vectors  such  that  the  vectors  fW\  ,  f W|  „  are  identical  and  Iwj  ,  , 

t  J  P  1  JP+2  l  JP+1 

{w}  are  identical  (the  subscript  p  designates  the  p^  iterative  cycle).  This 

case  occurs  when  the  matrix  being  iterated  upon  has  two  roots  equal  in  magnitude 
but  differing  in  sign.  The  relevant  eigenvalues  are  then  the  two  square  roots  of  the 
product  of  the  two  apparent  eigenvalues,  i.e., 


(A-7) 


The  above  condition  occurs  in  practice  when  the  applied  loading  is  a  bending 
moment  across  the  width  of  the  panel.  It  is  obvious  that  the  same  value  of  moment 
will  induce  buckling  independent  of  whether  the  moment  acts  in  either  of  two  oppos¬ 
ing  directions. 

The  analysis  procedure,  as  described  above,  yields  a  critical  stress  that 
corresponds  to  one  preselected  value  of  X  ,  the  buckling  wavelength  parameter. 
As  X  is  varied  the  critical  stress  will  vary,  reaching  a  minimum  for  an  as-yet 
unknown  value  of  X  .  To  establish  this  absolute  minimum,  one  can  first  select  a 
range  of  X  fs,  and  calculate  for  each  the  corresponding  critical  stress  parameter, 
n  .  Next,  employing  these  results,  a  polynomial  n  versus  X  relationship  is 
established  through  use  of  a  curve-fitting  technique.  The  minimum  n  is  then 
obtained  through  solution  of  the  equation  resulting  from  the  condition  d  fl 

d  X  “ 

The  value  of  the  critical  buckling  stress  parameter  is  a  function  of 


(1)  The  panel  stiffness  (represented  by  the  parameter  Dp). 

(2)  The  buckle  wavelength  (represented  by  the  parameter  X  ),  and 

(3)  The  applied  load  distribution  (described  by  the  parameter  p  ). 


Wavelength  considerations  are  taken  into  accound  by  the  sceme  described  above.  Thus , 
it  is  only  necessary  to  decide  which  ranges  of  the  parameters  Dp  and  p  are  of 

practical  interest.  A  discussion  of  this  question  is  given  in  Chapter  II. 


In  order  to  establish  effectively  convergent  results  all  computations  were 
performed  for  two  mesh  sizes,  the  first  being  the  9-interval  solution  described  here 
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and  the  second  an  18-interval  scheme.  In  general,  the  differences  between  the 
results  for  the  two  schemes  were  of  the  order  of  0.4%.  In  view  of  the  unavoidable 
inaccuracies  in  the  graphical  representations  which  have  been  made  of  the  computed 
results,  the  usefulness  of  any  improvement  of  the  results  through  an  extrapolation 
technique  is  negated. 
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APPENDIX  B 

PROCEDURE  FOR  DETERMINING  THE  STRESSES  AND  DISPLACEMENTS 
FORA  RECTANGULAR  SANDWICH  PANEL 


The  governing  small  deflection  theory  differential  equation  for  rectangular 
sandwich  panels  subjected  to  normal  pressure,  uniform  inplane  loading,  and  a  con¬ 
stant  temperature  gradient  through  its  thickness  has  been  developed  in  Reference 
22.  This  equation,  derived  from  elementary  equilibrium,  compatibility,  and  stress- 
strain  relationships,  is  applicable  to  a  sandwich  panel  having  unequal  shear  stiff¬ 
nesses  but  equal  bending  stiffnesses  with  respect  to  the  two  principal  directions. 
For  the  purpose  of  developing  the  design  curves  described  in  Chapter  III,  it  was 
deemed  desirable  to  consider  only  rectangular  sandwich  plates  having  isotropic 
cores  and  subjected  to  uniaxial  inplane  compression  and  uniform  normal  pressure. 
In  this  case,  the  differential  equation  reduces  to  the  form 


(1  -  /x)Ds 

77  D 
'  x  q 


,6  -,6  ,6  ,6 

,  o  w  .  d  w  /0  Jw  d  w 

<1+7x>TT  +  (1  +V  TT71 +  (3  V7T74  ~T 

0  x  ox  oy  dx  dy  dy 


Vx 


4  4  4 

(2  +  (3-  /x  )  77  )  d  w  +  (4  +  (3-  /x  )  77  )  d  w  n  +2  &  w 

x 


x  ax4 


2  2  7 

dx  dy  dy  -1 


2D.  ,2 

_ d _ w 

D  2 

s  dx 


2p 


D  17 
s  '  x 


(B-l) 


N 


where  'q  anc*  is  ^1C  inplane  loading  in  the  x  -  direction  (see  Figure  III- 1 

q 

in  Chapter  III),  p  is  the  uniform  normal  pressure. 


Since,  for  the  boundary  conditions  of  pi’esent  interest,  an  exact  solution  to 
(B-l)  cannot  be  obtained,  a  scries  solution  is  effected.  The  complementary  func¬ 
tion  portion  of  the  solution  ( w  )  is  taken  to  have  the  form 
CO  CO 

W  =  X  x  cos  - - —  +  X  Y 

C  b 

n=l  m=l 

where  X  is  a  function  of  x  only 
and  Y  is  a  function  of  y  only. 


cos 


m7rx 


(B— 2) 
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Substitution  of  Equation  (B-2)  into  (B-l)  yields  two  differential  equations. 


a-/oDs 


d6X 


x  q 


2(1-  )DS  /  n  7 r  \  ^  1 

(1+VV(  (;wx} 

x  q  < x 


n  d4X 


dx^ 


4* 


(1+n>t>s  (nr)  (3+7?  )  +  —  (— f(4  +  (3-M)’)  >+iik 

r)  D  \D/  x  *J7 _ '  k  '  x  D0 


x  q 


(1-  u  )D.  . 

^  s  /  n7T  \  + 


V  b 

4  "I 


i?  D 
L  x  q 


Vx  \  b 


n7T 


X  =  0 


d2X 


dx 


(B-3) 


and, 


(1- f±)D  d  Y  _ 


V  D  dy' 

'  x  q  J 


(1-  /a  )D  /mTT 


Tj  D 
'  x  q 


(3  +  77  )  +. 


x  V  x 


d4Y 

dy4 


s  /  m7r 


17  D 
x  q 


’  1  2 

)  (1  +  77  )+ - (JHZL)  (4+  (3-//.)  77 

/  x  77  x  \  a  '  x 


d2Y 


dy' 


6  4  2 

(1-  m)D  ,  m7r  \  ,,  4  1/  m7 t\  ,n  ,  2D  /  m7r 

r  s  '  — ]  I+77  +  —  f - ]  2+  3  -/i  77  )  +  _q( 

a  I  x  77  x\  a  /  v  D 


77  D 
x  q 


=  0 


(B— 4) 


Equations  (B-3)  and  (I>4)  ave  ordinary  linear  differential  equations  and  can 
be  readily  solved  to  yield  expressions  for  X  and  Y  .  Also,  the  particular  integral 
can  be  obtained,  using  procedui'es  described  in  Reference  22,  resulting  in  the 
following  complete  expression  for  the  deflection 
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+  2  K4  cosh  y  +  K5(<t>*-1)  cosh  4>5~-y 

m(odd) 


+  V*/-1'  “h*6-5rL 


a2  2 

- X 

4 


b2  2 
—  y 
4 


where  c/>  ^ ^  are  the  roots  of  the  cubic  auxiliary  equations  of  Equations 
(B-3)  and  (B-4)  and  K^, —  are  the  six  constants  required  to  satisfy  the  boundary 

conditions.  It  can  occur  for  some  combinations  of  dimensions  etc.  that  the  auxiliary 
equations  have  complex  roots.  In  such  cases  <£ ^  and  <p  ^  become  complex  conjugates 

and  is  purely  imaginary.  In  principle  the  hyperbolic  functions  with  complex 
0 

arguments  should  be  transformed  into  the  appropriate  trigonometric  and  hyperbolic 
functions.  However,  if  the  complete  computation  is  carried  out  using  a  complex 
mode,  as  is  possible  in  FORTRAN-coded  programs,  it  is  not  necessary  to  perform 
this  transformation  and  the  present  algebraic  form  can  be  retained  without  loss  of 
generality. 

For  the  special  case  of  =  0  (no  midplane  force),  equal  roots  occur  in  the 
two  auxiliary  equations,  so  that 

0!  =  tf>2  =  1  and^4  =  4>5  =  1 

and  the  deflection  function  has  the  slightly  different  form 


w  =  2 


n(odd) 


„  n7 rx  .  .  ii7rx  „  ,  n7rx  t„  ,  i  nvx  mry 

K — ; —  smh  — - +K  cosh—; —  +K  cosh<p„ — : —  cos - - 

01  b  b  01  b  03  3  b  b 


lllTT  y  .  ,  m  7T  X  ,  iuh  y  .  ,  m "  v  wo  in  " 

Kni - -  s“h_^r-  +  Knr  cosh  -IT-  '■  Knfi  c0sh  *  6— r1  — 


m7ry  cos  m7rx 


m  (odd) 

P  |  4  a4  4  b4 

h48Ds  [  (X  ~  16  '  ^  '  16  ) 
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The  solutions  presented  in  Chapter  III  can  be  divided  into  the  following  three 
classes  of  problems: 

(1)  Simple  support  on  all  edges  with  a  uniform  normal  pressure  and  a  constant 

temperature  difference  between  thefaces.  (17  =  0 .) 

(2)  As  in  case  (1)  but  with  the  addition  of  uniaxial  uniformly  distributed 

inplane  compression.  ( rj  ^  0). 

(3)  Simple  support  on  two  opposite  edges  and  fixed  support  on  the  other  two 
edges  and  with  same  loading  as  case  (1). 

The  boundary  conditions  pertinent  to  these  cases  are  as  follows: 

For  Cases  (1)  and  (2) 

at  x  =  ±  ~~  :w  =  0,M  =0,V  =  0 
2  x  y 

L. 

at  y  =  ±—  :  w  =  0  ,  M  =  0,.  V  =  0 
J  2  y  x 

where  V  and  V  are  vertical  shear  stress  resultants  in  the  core 
x  y 

For  Case  (3) 


at 

,  d. 

X  =  ±  — 

:  w  =  0,  M  =  0, 

M 

=  0 

2 

X 

xy 

at 

'< 

11 

14. 

>|cr 

:  w  =  0,  J3  =  0, 

M 

=  0 

y 

xy 

where  is  the  twisting  moment  and  /3  y  is  the  rotation  of  the  vertical 

fibers  of  the  core.  Note  that  although  the  classic  solution  for  an  isotropic 
plate  requires  only  two  boundary  conditions,  three  are  required  for  a  sandwich 
panel. 

Substitution  of  the  above  conditions  into  the  appropriate  expressions  for 
moment  and  shear  leads  to  simultaneous  equations  from  which  the  constants  (K^--^ 

or  K01 — KQg)  are  extracted.  The  deflections  and  moments  are  then  obtained  by 
back  substitution  into  the  pertinent  expressions. 

In  cases  (1)  and  (2)  the  series  involved  are  orthogonal  while  in  case  (3)  they 
are  non-orthogonal.  Thus,  in  cases  (1)  and  (2),  the  individual  constants  can  be 
determined  directly  for  any  given  term  of  the  series,  defined  by  m  or  n,  indepen¬ 
dently  of  any  other  term  of  the  series  and  the  problem  is  reduced  to  the  simple  solu¬ 
tion  of  two  groups  of  three  simultaneous  equations  in  Kp  K2,  and  K4,  Kg,  Kg. 

In  the  computer  program  employed  for  solutions  to  case  (1)  and  (2)  problems,  a  con¬ 
vergence  criterion  is  used  to  select  the  number  of  terms  in  the  series. 
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In  case  (3)  the  values  of  Kj — Kg  for  any  term  of  the  series  depend  on  all 
other  terms.  For  a  series  of  n  terms,  6n  simultaneous  equations  must  be  solved. 
Here,  since  the  computational  time  will  be  quite  long  for  a  series  with  many  terms, 
both  convergence  and  computational  efficiency  considerations  have  been  used  to 
govern  the  number  of  terms  selected. 
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APPENDIX  C 

SANDWICH  CYLINDER  INSTABILITY  UNDER  CIRCUMFERENTIALLY 
VARYING  AXIAL  STRESS 


This  problem  concerns  the  prediction  of  the  elastic  instability  of  a  circular 
cylinder  of  sandwich  construction  subjected  to  axial  stresses  that  vary  in  the  cir¬ 
cumferential  direction.  The  geometric  and  load  conditions  for  this  problem  have 
been  detailed  in  Chapter  IV  (see  Figure  IV- 1). 

The  solution  approach  is  essentially  the  same  as  that  employed  for  isotropic 
cylinder  analysis  in  Reference  9.  There,  Donnell's  partial  differential  equation  has 
reduced  to  an  ordinary  differential  equation  by  assuming  the  buckled  shape  to  be 
sinusoidal  in  the  axial  direction.  Then,  the  ordinary  differential  equation  was 
reduced  to  a  set  of  algebraic  equations  through  use  of  second-order  finite  differ¬ 
ences  and  matrix  iteration  was  applied  in  determination  of  the  eigenvalues  and 
eigenvectors. 


Now,  Donnell's  equation,  modified  to  apply  to  a  sandwich  cylinder,  becomes: 
(see  Reference  13). 


D 

s 


„  8  2 
V  w  +  2tr 


Ds  2 


i  -  ttS  V 
D  r^ 

q 


(e  i-S  + 


4 

v  cr 

x 


=  0  (C-l) 


where  x  in  this  case  is  a  non-dimensionalized  axial  coordinate  (the  axial  coordinate 
divided  by  the  radius). 


The  displacement  w  is  assumed  to  be  of  the  form 


w  =  W  sin  X  x 


(C— 2) 


where  W  is  a  function  of  the  circumferential  coordinate  (y)  only,  and  X  is  a 
function  of  the  wavelength  in  the  axial  direction. 


Substituting  Equation  (C-2)  into  (C-l)  and  introducing  a  nondimensionalization 
of  the  stress  distribution  similar  to  that  used  for  the  flat  panel  in  Chapter  II  and 
Appendix  A,  the  differential  equation  reduces  to 
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T  fi  24  9  ^  X  6  9  ^  X  IT 

(Dg  X  +  2  tr  E  X  +  2  tr  — —  )  W  -  (4  Dg  X  +  2  tr  — - — )  W 


„  „  v  4  ...IV  ,  _  ,2  ...VI  _  ...vm 

+  6  D  X  W  —  4  D  X  W  +DW 
s  s  s 


-  2  tr  cr 


6  \  2  4  9  \  2  IT 

.  X  (1+— )  (^SW)-  X  (2+—  )(PW)1 


+  .2(1  +  3^)( 

C  c 


(  Ps W)VI  ]  =  o  (C-3) 


where  primes  indicate  derivatives  with  respect  to  y  only  and  = 


The  conditions  examined  are  limited  to  stress  variations  which  are  symmet¬ 
ric  about  the  vertical  (z)  axis.  In  order  to  write  Equation  C-3  in  terms  of  finite 
differences,  the  circumference  of  the  cylinder  is  divided  into  2n  spacings,  each  of 
length  rT^n  .  as  shown  in  Figure  IV-1-  Note  that  elements  0  and  n  lie  on  the  axis  of 

symmetry.  With  use  of  second  order  finite  differences,  C-2  reduces  to 


Wk-3  +  %3)  AI 


II  (  ^k-4Wk-4  +  ^k+4Wk+4  B1  (  ^  k-3Wk-3  +  ^k+3Wk+3 )  B2 

+  (Pk-2Wk-2+^k+2Wk+2  )  B3  -(^k-lVl  +^k+lWk+l)B4 

+  ^,WkB6]  =0  (C-4) 


-  /  w,  +  w,  \ 

A  + 

I  W,  W,  / 

\  A 

\  k-5  k+5  / 

1 

1  k-4  +  k+4 j 

1  2 

+  ( w,  +W,  ) 

A,  - 

(  W,  +W,  ) 

A 

'  k-2  k+2' 

4 

'  k-1  k+1 ' 

5 

where 


V^T  A2=a3(T^+1> 

A  =  A  2 (23  A  2  +  12A  +  1) 

u 

9 

3  2  1  7T  KA 

V  a  (C8A  4  52  A  4  12  A  +y  > 

A  =  A  (126  A  3  +  116  A  2  +  39  A  +IL  +  L ..L_KA) 

5  3  3n2  D 
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APPENDIX  D 

THERMAL  STRESSES  AND  BENDING  MOMENTS  IN  A  HEATED 
SANDWICH  CYLINDER  SUPPORTED  BY  RIGID  BULKHEADS 


A  sandwich  cylinder  of  cross-sectional  properties  identical  to  those  assigned 
to  the  sandwich  cylinder  of  Appendix  C,  but  of  finite  length  2c,  is  supported  at  each 
end  by  a  rigid  bulkhead.  The  cylinder  may  be  continuous  over  several  bays,  so  that 
each  bay  is  effectively  built-in  at  the  rigid  bulkheads  due  to  symmetry,  or  one  bay 
alone  may  be  considered  with  simply  supported  ends.  The  cylinder  is  uniformly 
heated  so  that  a  temperature  change  from  the  stress-free  state  (A  T)  is  sustained, 
while  the  bulkheads  are  unheated. 


The  subject  problem  is  similar  to  that  of  a  beam  on  an  elastic  foundation. 
Using  this  concept,  the  basic  differential  equation  for  the  radial  deflection  (w)  is 


d4w  2tE  d2w 

d^  r2D  dx2 

Q 


+  —  (w  +  r  a  A  T)  =  0 
r2Dg 


(D-l) 


where  a  is  the  coefficient  of  thermal  expansion.  A  general  solution  for  this 
equation  is 


w  =  C^  cosh  m 


Et 
Hr2D 


1/2 


x  +  C2  cosh  m2 


Et 

r2D 


1/2 


x  -  raAT  (D-2) 


where  m^  and  the  roots  of  the  auxiliary  equation,  are  given  (in  nondimensional 
form)  by 


nr 


with 


1  +  *J\ 


EtDs 

r2  D  2 

q 


and,  mr 


d  J 
Et 


1  -  y/l 


(D-3) 


DcDq 


The  roots  and  m2  may  be  real  (H^  >  2)  or  complex  (H^  <  2).  In  the 
latter  case,  the  hyperbolic  functions  of  Equation  (D-2)  should  be  modified  as  follows: 

cosh  (a  +  i  b)  =  cosh  a  cosh  b  +  i  sinh  a  sin  b 

but,  as  in  Appendix  B,  provided  all  computations  are  carried  out  in  a  complex  mode, 
it  is  not  necessary  to  transform  the  hyperbolic  functions. 
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As  noted  above,  two  support  conditions  -  simple  and  fixed  support-- are  of 
present  interest.  Consider  first  the  simple  support  condition,  where  the  boundary 
conditions  are 

d^w  2tE 

at  x  =  ±  c,  w  =  0  and  M  =  -D - (w  +  r  a  A  T)  =0 

x  s  2  2™ 

dx  rzD 


Using  these  conditions,  the  constants  and  C2  are  determined  and  the  deflected 

form  becomes  _  ,  , 

^  Am  „  cosh  m  H  „  cosh  m  H  „  n 

r  C£  A  T  ,  2  lx  ,  2  2  x  ,  2  2, 

w  = - (mn  -2  — : - - - (nv  -2  - ; - —  -(m„  -m,  ) 

n  n  'O'  OapK  vn  U  *  1  '  oaoU  m  14  '  O  1  ' 


(m2  ) 


2  2,  2  cosh  m  H  1  cosh  m  H  2  1 


where  H  =  ~  —  and  H  =  — 

x  \Dq/  r  1  \Dq/  r 


(D-4a) 


The  formula  for  hoop  stress  is 


E  T  *  _  "1 

=  -  -  w  +  raAT 
r  J 


Thus,  the  hoop  stress  at  the  center  (at  x  =  0)  is  given  by 


E.a  AT 


m  2-2 


m^-2 


<£  ,2  2,  cosh  m  H.  cosh  m  H. 

r  c  (m  -m  )  L  11  2  1 


(D-5a) 


Also,  the  expression  for  longitudinal  bending  moment  is 


M  =  -  D  —  -  —  (w  +  raAT) 
X  S  dx2  r2D_ 


Hence,  the  maximum  bending  moment  in  the  cylinder  (at  the  center,  x  =  0)  is 

On  *-!)(»,*-!)  f  j  i  ", 


=  H  .  D  r  a  A  T 


(mi2-m22) 


cosh  m  H  cosh  m  11 

.LI  ^  I 


(D-  6a) 


With  built-in  ends,  the  boundary  conditions  arc 
,  dw  Qx 

at  x  =  1  c,  w  =  0  and  —  = 

dx  Dq 


where  Q  is  the  transverse  shear  force  and  is  given  by 


Q  =  -D 

x  s 


d  w  2tE  dw 

dx3  „2n  dx 
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APPENDIX  E 

BUCKLING  OF  A  SANDWICH  CYLINDER  UNDER  THERMALLY 
INDUCED  HOOP  STRESS 


1,  GOVERNING  EQUATIONS 


Appendix  D  treats  the  problem  of  determining  the  hoop  stresses  in  a  uniformly 
heated  sandwich  cylinder  supported  at  either  end  by  rigid  unheated  bulkheads.  This 
Appendix  develops  a  method  for  determining  the  cylinder  temperature  change,  Atci>, 
for  which  the  cylinder  will  buckle  in  consequence  of  the  induced  hoop  thermal 
stresses,  Thus,  it  is  assumed  that  the  hoop  thermal  stresses  will  have  already  been 
determined  by  use  of  the  technique  of  Appendix  D. 


For  the  present  case,  Donnell's  equation,  modified  to  apply  to  sandwich  con¬ 
struction,  takes  the  form  (Reference  13). 


D  V  w  +  2t 
s 


r  dA  a  A 2 

E  w  „  4  0 

—  - I+V  <r  — 2 

r  dx  y  dy 


=  0 


(E-l) 


form 


The  general  solution  for  hoop  stress,  <jy,  as  derived  in  Appendix  D,  is  of  the 


cr  =  Acoshm,H  +  B  cosh  m„H  (E-2) 

y  lx  2  x  ' 

As  described  in  Chapter  VI,  three  separate  cases  must  be  considered: 

(a)  Fixed  support  condition.  In  this  case  the  cylinder  consists  only  of  one  bay 
fully  fixed  against  rotation  at  each  rigid  bulkhead. 

(b)  Simple  support  condition  -  one  bay.  Again  the  cylinder  consists  only  of 
one  bay  but  the  ends  are  no  longer  restrained  against  rotation  at  the 
bulkheads. 

(c)  Simple  support  condition  -  many  bays.  The  cylinder  is  taken  to  be  con¬ 
tinuous  over  many  bays  as  in  a  fuselage.  The  bulkheads  provide  no 
restraint  against  rotation  but  restrain  completely  all  radial  displacement. 

In  Chapter  VI  it  was  demonstrated  that  case  (c)  is  derived  from  a  combination 
of  fixed  support  thermal  stresses  and  a  simply  supported  buckling  mode.  Thus  no 
theoretical  treatment  of  this  case  is  necessary  in  this  Appendix  and  only  cases  (a) 
and  (b)  are  discussed  here. 
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2.  FIXED  SUPPORT  CONDITION 

To  analyze  for  fixed  support  conditions,  the  origin  of  coordinates  is  first 
shifted  from  the  center  of  the  cylinder  to  the  left  bulkhead  (see  sketch  below).  The 
hoop  stress  equation  then  takes  the  form 

0‘y  =  Ea  At[  A  cosh  m^  (Hx-Hj  )  +  B  cosh  m2(Hx-H^ )]  (E-2a) 


Also,  since  the  hyperbolic  functions  are  not  convenient  for  use  in  Equation  (E-l),  a 
transformation  to  an  infinite  cosine  series  is  accomplished  with  Fourier  analysis 
techniques,  resulting  in 


c  =  E  a  A  T  I  P 

n  1  +  8, 
p=0  Op 


where  8  is  the  Kronecker  delta  /  8  \  \  =0,  i  /  j 
Op  1J 

\8ij  =  1'  i=J 


and  C  = 
P 


2  Am^Hj  sin  m-^H^ 


—  )  +  (m^j) 


2B  n^Hj  sinh  m2^l 


(m2Hl 


(E— 3) 


(E-4) 


for  p  even,  and 


To  achieve  a  solution  for  Equation  (E-l)  with  use  of  Equation  (E-3),  Galerkin's 
method  is  employed.  Thus,  a  deflection  function  which  satisfies  both  the  static  and 
kinematic  boundary  conditions  is  selected.  This  function,  which  was  employed  by 
Anderson  in  Reference  28  is 
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w  =  sin 


7T  y  V  ,  1V  7TX  .  7 rx 

IT  ^  am  cos(m-l)— -  cos(m+l)  — 

A  m=l  L 


7r  r 

where  X  =  -^-is  the  half  wavelength  of  the  buckled  shape  in  the  circumferential 

direction  and  the  Coefficients  am  range  from  m=l  through  m=k  (the  function  for  w  is3 
of  course, an  infinite  series  with  k  =  CO  but  is  truncated  at  a  finite  value  of  k  for 
numerical  analysis).  The  unknown  coefficients  am  are  determined  from  the 
condition 


X  „  2c 


w  •  Q  (w)  dxdy  =  0 


where  Q  (w)  is  the  differential  equation  (E-l)  with  cr^  from  Equation  (E-3). 

After  integration  of  Equation  (E-6)  and  rearrangement  of  terms  one  obtains 

fm[V1  +  Slm>  + V2]-V2Um  <X-  S2m>-am+2  V2 


+  ^  ^  ai  \Cm-i-2  '  2Cm-l  +  Cm-i+2'  °m+i-2  +  2  Cm+1  ’  ^+2^2 


Hd  ,f[(l+  /52)(m-i 


-i-2)2  +  2  (i+ 1) (m-i-2)  +  (i+1)2]  Cm_._ 


-2  [(1  +  /3 2)(m-i)2  +  2i(m-i)  +  (i2  +  1)  ]  Cm_. 

+  [(1  +  /3  2)(m-i+2)2  +  2(i-l)(m-i+2)  +  (i-1)2  j  Cm_.+2 

-  [(1+  /32)(m+i-2)2-2(i-l)(m+i-2)+  (i-1)2]  Cm+._2 
+2  [(1  +  P2)(m+i)2  -  2i(m+l)  +  (i2  +  1)  ]  Cm+1 

-  [(l+/92)(m+i+2)2-2(i+l)(m+i+2)+(i+l)2]  Cm+i+2 
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where  U 

m 


r  2  2 1 

[(m-1)  +/8 


2  +34^4 .  C”-1)4  l2 ,8H12  (m-1)4 

T  Hd  [  (m-1)2  +£2]  T  2  [(m-l)2+^2j 

and  n  =  °  1  C°cr 


Equation  (E-7)  can  now  be  written  for  m  =  1,  2,  3,  —  k,  which  yields  k 
simultaneous  equations.  The  condition  of  a  vanishing  determinant  leads  to  the  criti¬ 
cal  value  of  the  buckling  stress  parameter  17  .  It  is,  however,  more  convenient  as 
in  Appendix  A,  to  present  the  problem  in  matrix  form  and  apply  an  iterative  tech¬ 
nique  for  the  determination  of  FI  for  a  number  of  values  of  the  wavelength  para¬ 
meter  /3  .  When  the  simultaneous  equations  are  written  in  matrix  form,  one  obtains 


[u]  {a}  ♦  n  [I]  {a}  =0 

ifM  =  ■  [U1  l1  1  { a )  (E'8) 

It  is  found  that  the  matrix  [  U  ]  *  [  £  ]  has  a  checkerboard  pattern,  i.e., 

alternate  terms  are  zero.  This  matrix  can  be  rearranged  to  form  two  decoupled 
systems,  corresponding  to  symmetric  and  antisymmetric  failure  modes,  which  can 
be  treated  separately.  The  symmetric  mode  leads,  in  all  cases,  to  the  lower  criti¬ 
cal  stress.  From  the  eigen  value  IT  of  the  system  the  value  of  Cocr  =  C0  E  a  A  Tcr 
can  be  computed  leading  to  the  critical  temperature  rise  Al^. 


Unfortunately,  this  approach  demonstrated  such  poor  convergence  that  it  was 
not  possible  to  determine  the  critical  temperature  even  with  m  =  30.  As  noted  in 
Chapter  VI,  however,  a  more  critical  mode  of  buckling  than  the  present  case  can  be 
defined  (Case  (c))  when  the  intent  is  to  analyze  a  long  cylinder  supported  over  many 
rings. 


3.  SIMPLY  SUPPORTED  EDGES  -  ONE  BAY 


As  in  the  case  of  built-in  edges  the  expression  for  the  hoop  stress  in  the 
cylinder  is  transformed  into  the  Fourier  cosine  series  (Eq.  E-3),  where  the  coef¬ 
ficients  Cp  are  given  by  Eq.  E-4. 

The  analysis  follows  that  of  Hoff  (Reference  14),  suitably  modified  to  include 
shear  deformation  terms.  Hence,  the  displacement  function  is  taken  to  have  the 
form 


„  mrr  x 

w  =  cos  ny  2.  a  cos  - - — 

m  \ 


(E— 9) 
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Substitution  of  Equations  (E-3)  and  (E-9)  into  (E-l)  yields  an  expression  in¬ 
volving  double  series.  These  series  can  be  manipulated,  as  in  Reference  14,  to 
generate  a  set  of  linear  simultaneous  equations  (see  Equations  28  of  Reference  14), 

The  general  form  of  one  such  equation  is  the  counterpart  of  Equation  (E-7)  of  the 
fixed  support  case  and  the  matrix  form  of  these  equations  has  the  same  appearance 
as  Equation  (E-8).  The  expression  for  Um  given  by  Hoff,  however,  must  now  include 
the  shear  deformation  terms,  i.e. 


In  the  present  case  the  matrix  associated  with  this  term  (  [u]  )  is  a  diagonal 
matrix  and  is  therefore  easily  inverted.  Using  the  matrix  iterative  techniques  of 
Appendix  A  the  eigenvalues  of  the  system  for  a  range  of  values  of  the  wavelength 
parameter  are  computed  and  the  minimum  determined  using  a  standard  curve 
fitting  procedure. 

The  eigenvalue  thus  determined,  R  in  Hoff's  notation,  corresponds  to  the 
critical  thermal  strain  d  A  Tcr. 
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APPENDIX  F 

SANDWICH  CYLINDER  STRESSES  DUE  TO  INTERNAL  AND  EXTERNAL 
PRESSURE  AND  RADIAL  TEMPERATURE  GRADIENTS 


The  cylinder  treated  in  this  section  is  of  the  identical  cross-section  as  those 
examined  in  Appendixes  C,  D,  and  E  and  is  infinitely  long.  Now,  however,  a  dis¬ 
tinction  must  be  made  between  the  radii  of  the  inner  and  outer  skins.  The  inner 
radius  is  designated  as  r^  and  the  outer  radius  as  rQ.  The  inner  and  outer  skins 
are  at  constant  temperatures  T^  and  T0,  respectively  (AT  =  T0-Tj),  and  both 
internal  and  external  pressures  may  be  present.  Also,  the  core  is  now  assumed 
to  have  a  modulus  of  elasticity  in  the  radial  direction  (Er)  and  zero  stiffness  in  the 
circumferential  and  axial  directions. 

In  this  approach  to  the  development  of  relationships  for  the  stresses  in  the 
faces  and  core,  each  cylindrical  component  (the  inner  and  outer  faces  and  the  core) 
is  examined  separately  and  finally  combined  using  conditions  of  compatibility  at  the 
interfaces.  Axial  stresses  are  transmitted  from  one  cylinder  to  the  next  only  at 
the  ends  of  the  cylinders;  such  "end  effect"  stress  systems  are  not  considered  in 
this  development. 


Sketched  below  is  a  representation  of  the  forces  acting  on  the  skins  and  on 
the  core.  The  unknown  interface  stress  between  the  inner  skin  and  the  core  is 
designated  asqj .  At  any  radius  r,  the  radial  stress  within  the  core  must  be 

q.r. 

q=  -7-  (F-D 

ri  qi 

and  it  follows  that  q0  =  Qi  “  =  —  (F-la) 

r  0 

where  R  =  ,  (F-2) 


Inner  Face 


Core 


Outer  Face 
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while  the  radial  strain  due  to  this  stress  is 

q.r. 

*  c  -  17  <M> 

By  integration  from  rj  to  r0,  one  obtains  the  radial  growth  (i.e.,  the  change  in 
radial  displacement  between  the  inner  and  outer  faces)  due  to  this  strain  as  follows: 

A  w  =  /  =  q  ri_  ln(R)  (F— 4) 

4  r.  E  r  E 

1  r  r 


The  radial  growth  due  to  thermal  strains  is  simply 
«...  _/To+Ti\  ,  . 


-PP) 


(ro“ri) 


(F-4a) 


Considering  now  the  total  radial  growth  of  the  core  (wc)  as  being  composed 
of  both  the  growth  due  to  thermal  expansion  and  A  wc  ,  one  obtains 

A  wc  =  [-2(VTi)<Vri)  +  I7  ln(R)]  ,F"5 

Radial  displacement  of  the  inner  cylinder  (Wj)  is  given  by 


'i  =  rj  [a  Ti  + 


^1  +  Pi)  ri  crx. 

-  “  r1  — - 


where  17 x.  is  the  axial  stress  in  the  inner  cylinder,  while  the  strain  in  the  axial 
direction  is 

r  (Qi+Pilri  I 

=  la  Ti  +  T  -/x  J  (F~ 


Similarly,  for  the  outer  cylinder,  the  radial  displacement  and  axial  strain  are 

r  „  „  %  (%>Poiroi 


r  xo 

r  =  r  I  a  T  -  H-  — —  - 
0  0  L  0  E 


=  fa  t  +  +  n  (qo+Po)rol 

L  0  V  Pf  J 


where  cr  the  the  axial  stress  on  the  outer  cylinder. 
xo 
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